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Abstract 

It is well-known that there is a geometric correspondence between high-frequency quasi-normal 
modes (QNMs) and null geodesics (spherical photon orbits). In this paper, we generalize such 
correspondence to charged scalar field in Kerr-Newman space-time. In our case, the particle and 
black hole are all charged, so one should consider non-geodesic orbits. Using the WKB approxi¬ 
mation, we find that the real part of quasi-normal frequency corresponds to the orbits frequency, 
the imaginary part of the frequency corresponds to the Lyapunov exponent of these orbits and 
the eigenvalue of angular equation corresponds to carter constant. From the properties of the 
imaginary part of quasi-normal frequency of charged massless scalar field, we can still find that the 
QNMs of charged massless scalar field possess the zero damping modes in extreme Kerr-Newman 
spacetime under certain condition which has been fixed in this paper. 
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I. INTRODUCTION 


In 1970, Vishveshwara found that there exist complex frequency modes in perturbations 
of Schwarzschild black holes[T], which are called quasi-normal modes (QNMs). QNMs play 
an important role in the physics of black holes [2], where many works have been done[3HH]. 
Various aspects of QNMs continue to be uncovered in recent years. Actually, there have 
been some nice reviews of this subject, describing the methods to calculate QNMs and in¬ 
troducing the developments of QNMs: Ref. [9] shows some properties of eigenfrequencies and 
potential applications of QNMs to gravitational wave asteroseismology; Ref. HDl introduces 
some analysis about the application of QNMs to relativistic stars and the detection of QNMs; 
Ref. [H] discusses the applications of QNMs in the AdS/CFT duality and to astrophysical 
black holes; Ref. na discusses the gravitational instabilities in higher than four dimensions 
and the AdS/CFT interpretation of QNMs. 

The semi-classical interpretation of QNMs is particularly interesting. In 1984, Ferrari 
and Mashhood[T3] found that the quasi-normal frequency in Schwarzschild black hole can 
be written as 


u ^ 





IL 


when Z 3> 1. Here H is the frequency of circular photon orbits and 7 l is the Lyapunov 
exponent of the orbits which describe how quickly the cross section of the null geodesic con¬ 
gruence changed under inhnitesimal radial perturbations. In the eikonal limit, the authors 
of IH] studied the QNMs by the WKB approximation and obtained the expression of the 
imaginary part of the quasi-normal frequency. The relationship between the QNMs and the 
circular null geodesics was then studied in [15] , which showed that the QNMs of black holes 
in any dimensions are determined by the parameters of the circular null geodesics. Recently, 
the authors of [T6| made a detailed study of the geometric interpretation of the QNMs of 
massless scalar field in general Kerr black holes by the WKB approximation. 

Usually, people may think that the charge effect is unimportant since most celestial bodies 
are electroneutral. But things are not always like that. Recently, some super-Chandrasekhar 
white dwarfs have been found [191 EO] and this phenomenon can be explained by the charge 
effect [2TII23] . Even if the particle is uncharged, the behaviors of its trajectories near the 
singularity of Schwarzschild and Reissner-Nordstrom (RN) black holes are quite different. 
For the RN black hole, the singularity is untouchable by geodesics, but this phenomenon 
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does not occur for the Schwarzschild black hole. Furthermore, when the matter helds or 
perturbations in Kerr-Newman (KN) black hole are charged, the interaction between matter 
helds and electromagnetic held in the KN black hole will have signihcant inhuence on the 
QNMs, which has been studied in, e.g., [231 - I28] . 

So it is interesting to study the geometric interpretation of QNMs of the charged scalar 
held in the KN black hole, which is the main motivation of this paper. In order to achieve this 
geometric interpretation, we also need the behavior of particle motion in the KN spacetime, 
which can be found in [2SH3I], and generalize the method in (TU] to explore the geometric 
correspondence between the QNMs and spherical orbits of charged particles for the KN 
black hole, where special attention is paid to the non-geodesic nature of these orbits. It is 
shown that, noting the U{1) gauge invariance of physical quantities, the QNMs have perfect 
geometric interpretation by means of the non-geodesic world line congruence. The influence 
of charge of the particle is also addressed. 


We organize the rest of the paper as follows. In Sec. Ill we will investigate the Hamilton- 
Jacobi formalism of particle motion and study the geometric correspondence of QNMs. We 
can verify that those relationships [I6] also hold: in the leading order, the real part of quasi¬ 
normal frequency corresponds to the energy E, the azimuthal quantum number corresponds 
to the angular momentum Lz, and the real part of angular eigenvalue corresponds to the 
Carter constant Q; in next-to-leading order, the imaginary part of quasi-normal frequency 
corresponds to the Lyapunov exponent of radial motion, and the next-to-leading order 
correction of the angular eigenvalue corresponds to the imaginary part of the Carter constant. 
In Sec. |IV| we obtain the geometric correspondence by some specihc example and study the 
influence of charge. We also study and find the condition of the zero damping modes 
in extreme KN black hole. We can observe that the imaginary part of the quasi-normal 
frequency will approach to zero under this condition. In Sec. |V] we conclude with some 
discussion. 


5 



II. QUASI-NORMAL MODES IN KERR-NEWMAN SPACETIME 


We consider the charged massive scalar held u in KN spacetime in Boyer-Lindqnist co¬ 
ordinate system: 

2Mr-Q^' 


ds‘^ = 


P" 


dt"^ + ^dr'^ + p'^dd'^ + ( 1 ) 


2a 


Aa^ sin^ 9 j sin^ Odcf^ - - (2Mr — sin^ Odidcf). 


P 

with M is the mass of black hole, aM is the angnlar momentnm of black hole, Q is the 
black hole charge, cos^ 0, A = — 2Mr + and the fonr potential 

= —^[{dt)a — asiv? 9{d4))o\. The eqnation of motion of u is the Klein-Gordon(KG) 
eqnation with charge in cnrved spacetime: 


P (Va - iqAa){Vb - iqAb)u - = 0. 


( 2 ) 


with q is the charge of the scalar held, p* is the mass of the held. We can obtain that 
VaA“ = 0, A°'Aa = — A^ = + a^), A^ = ^ nnder this gange. We can separate 


A 


A 


the scalar held as: 


u[ 


l,m 

and obtain the radial eqnation and the angnlar eqnation as: 
d f dR\ ( {K + qQr)'^ 


(3) 


dr \ dr ) ^ 


A 


— [Aim + — 2.amuj + p*r^) i? = 0; 


1 d 


sinO^Uf)^ + — a^p^) cos^ 6* 


m 


+ Aim ) Mg — 0 


(4) 


(5) 


sin 9 d9 d9 J ' sin^ 6* 

with K = —uj{r^ + a^) -f- am, Aim is fhe angnlar eigenvalne of the angnlar eqnation and 
depend on qnantnm nnmber I and m. This is the Tenkolsky eqnation([32]) with spin s and 


sonrce are zero. 


A. The angular eigenvalue equation 


Under the transformation 


a; = In ( tan - ) , 


( 6 ) 
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namely dx = the angular eigenvalue equation O can be written as 


d'^ue 

dx'^ 


+ V^ue = 0 . 


( 7 ) 


with — fil) cos^ ^sin^ 6 — rn? + Aim sin^ 6. We can see that the influence of mass 

can be reflected in — /i^. We assume that m 7 ^ 0 and the condition 0 < /r* < |a;| holds in 
the following discussion. Since the quasi-normal frequency is complex with u = ur — iui, 
the angular value Aim is also complex: 

Aim = Aj^ + iAjm- ( 8 ) 

We can use the WKB analysis to deal with this equation, the approximate solution is 
c_ exp ^J-V^ {x') dx'^ 


Uq — 


V\x') 


X < x_ 


ue = 


U0 = 


a+ exp (^i \/V^ (x') dx'^j + a_ exp \/V^ (x') dx' 

I V<^{X') I 3 

c+ exp (^- ^J-V^ (a;') dx^ 

I V<^{X') I 3 ’ 


X > X^ 


(9) 

X- < X < Xj^ (10) 

( 11 ) 


where x_ and x^ is the turning point of V^{x), i.e. V^{x) is equal to zero at these two 
points. The wave propagates in the region X- < x < x+ and decay to zero when x —)■ ± 00 . 

The Bohr-Sommerfeld condition is 

rd+ 




de y a 2 {uj\ - cos 2 9 ^.^2 ^ 


+ ^fm = iL-\m\)Tr. 


( 12 ) 


with L = / -|- 1. We define fi = m/L, aR^a,^) = Af^/L"^, m* = Qr^q,^) = ojr/L. 

We can treat — ml as a small parameter and expand this integration to obtain 

9-1) 


CtR 


1 - 


2 L 2 


-(l-/r 2 ). 


(13) 


Using the facts that ojr ~ 0{l), uj ~ 0(1), m ~ 0{l) (see [IE]), we obtain that ur ^ U[. 
Then we can use the perturbation theory of eigenvalue equation [16], which leads to the 
next-to-leading-order correction of Aim- 


Aim = —2a^URUi{cos^ 9) = acui 


dAfm(z) 

dz 


z=aujji 


Finally we can obtain the approximate expression of Aim in eikonal limit {I 3> 1) as 


A 


Im 


/(/ + i) 


-h*) 


m 




(14) 


(15) 
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B. The radial eigenvalue equation 


We define 


Ur = A/r2 + a'^R, ^ 

dr^, dr 


and 
1/^ = 


A 


(r^ + 


^Im + fJ'-k'f' + 


2 2 A + 2r (r — M) 


3Ar^ 


+ a; 


ma + qQr 

r2 _|_ q 2 


r2 + a2 (r2 + a2)2j 

with Aim = “ 2amuj. Then the radial eq.@ can be transformed as 

d^u 


( 16 ) 


(17) 


drl 


+ V^Ur = 0. 


(18) 


In the condition of eikonal limit (Z S> 1), we can neglect the terms 
the expression of the effective potential V'', we rewrite the effective potential as 
[a;(r^ + a^) — ma — qQr^ — A{Aim + — 2maijj + 


m 


V’’ = 


(19) 


(r^ + a^)^ 

Under this transformation, when r —>■ +cx), we have r* —)■ +cx) and the effective potential 
—)■ — fil; when r —)■ r+, we have r* —>■ —cx) and —>■ (a; — )^- In order to 

obtain QNMs, we require that the wave is outgoing at r* —)■ +cx) and ingoing at r* —>■ —cx) . 
By comparing handling with the stationary Schrodinger equation, we know that boundary 
condition require that there is a point such that the effective potential is zero here. 
We can expand effective potential W here as 


V'(r,) = V„ + r„'(r. - rS) + ^(r. - rff + 0((r. - 

92 y 


( 20 ) 


with Uq = 0, Vo = W(rQ), Vq = Taking this Taylor expansion into 

radial equation and using the result of Schutz and Will[Tl] as well as boundary conditions, 
we can obtain 

1 


n + - = -: 

2 (2V'„")! 


( 21 ) 


Due to the spirit of QNMs, we write u: = ur — iuj, where the imaginary part is a small 


parameter in contrast to the real part. When we take this into (21) and Vq = 0, and expand 
with respect to tn, we obtain 

aw , 


W(r*,a;^)=0 = 




( 22 ) 


oJi = {n + -) 


1, dri 


2 ' ('dUi'i * 

1 di^ )rQ,0JR 


(23) 





















In other words, there is a radius ro such that 


dV^ 

I/^(ro,a;R) = 0 = 


(24) 


1 )^0,' 
= (^+ o)^aK- 


2^ m 


(25) 


duj Jro,ujR 


Then we can use (24) to calculate the position of apex of effective potential and the real 


part of frequency, and obtain the imaginary part of frequency. 


III. QNMS AND THE MOTION OF PARTICLE IN KN SPACETIME 


A. Geometric optics in KN spacetime 


The massive charged scalar held in curved spacetime satisfy the KG equation 


5'“^(Va - iqAa){Vb - iqAb)u - niu = 0. 


We assume that the wave function can be expressed as 


u = Ae^^. 


(26) 


(27) 


When we taking this formula (27) into the KG equation.(26) and dehne the 

leading order equation is 


g'‘“(k - qA ){R - qA^) + 


(28) 


and the next-to-leading order equation is 


2(F - qA'^)dfq In A + V^(/c'^ - qA>^) = 0. 


(29) 


We know that the mechanical moment kn — qAn is how wave propagate in curved background 


from leading order 28 and we note = k^ — qA^ to represents the four-velocity. The 


electromagnetic held makes particle’s worldline no longer geodesic and that the phase will 
change along the particle’s worldline: 


{U^ - qA^^)d^^ = -qii + 


p2A 


(30) 
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Furthermore, we can obtain U^'W from the leading order, that is the 
equation of motion of a charged massive particle in a curved spacetime with a electromagnetic 
held. It is interesting to know that the mechanical moment is different from the canonical 
momentum k^. As we all know that the electromagnetic held is gauge held and possesses 
the U{1) gauge invariance, when we do the gauge transformation 

^ + dfj.X- (31) 

then the phase of wave function will changed correspondingly: 

^ + X, (32) 


and the canonical momentum will change like the following: 


k'^ = ki, + d^x- 


(33) 


We can obtain that the wave front will change under gauge transformation but the four 
velocity will invariant = U'^. 


The next-to-leading order equation (29) gives rise to 


= -2U>^dfx\nA. 


Similar to [16], this equation is related to the huid expansion equation 




(34) 


(36) 


with ^ represents the area of cross section of the congruence. So the right hand side of 
above equation shows the change rate of cross section along the worldline. After comparing 
these two equations we have 


U>^d^\n{£/-2A) = 0 , 


(36) 


SO 

Aoc^/“T (37) 

It is well-known that the square of the length of amplitude represents the probability of 
the appearance of particle, this result tells us that the more the area of cross section of 
congruence, the less the probability of the appearance of particle. 
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B. Equation of motion for particles in KN spacetime 


We will use the Hamilton-Jacobi formalism to study the motion of charged massive par¬ 
ticle: 

- qA^){d,S - qA,) + = 0. (38) 

with S represents the Hamilton principal function. With two Killing vector fields and 
S can be written into the following form: 

S=^-Et + L,<t) + Sr{r) + Se{e). (39) 

where E represents the total energy of particle with gravitational as well as electromagnetic 
energy, Lz represents the z-directed specihc angular momentum. The property of timelike 
and axial Killing vectors ensure us that E and Lz are conserved constants. And r is the 
proper time of particle. The angular and radial equations are 

(S';)2 = Q- cos^ e (- a^E^ + , (40) 

\sm 0 J 

= {E (r^ + a") - aLz - qQrf - A ((L, - aEf + Q + , (41) 


where S'q represents derivative with respect to 9 and S'(. with respect to r, Q is Carter constant 
which is the third conserved quantity. We can see how the charge and mass influence the 
radial equation {qQr and fiY), but there is no charge influence in angular equation. We 
dehne 


7^(r) = {E Y + «') - aLz - qQr^ - A {{Lz - aEf + Q + YY) , 


Then we have 


0 ( 0 ) = Q-cos^e 


Sr(r) = 


r2 

^2 172 I „ 2,,2 

. 2 a ~ a E + a 
sm 9 

r Vn 


A 


dr 


Sg{9) = j y/ed9. 

So the Hamilton principal function has the form 

Vn 


s = 




Et J- Lz(j) T 


A 


dr + y/O d9. 


From ^ = 0 and |^ = 0, we can obtain 


rr ^2 


rd ^2 


y/n 


dr + 


cos^ 9 


d9 = r, 


(42) 

(43) 

(44) 


(45) 


(46) 
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: dr = 


-.d9. 


Vn'^' J Vq 

Taking derivatives of these two equations with respect to A: 


we have 


dX ^ dr’’ 


dr a? cos^ 9 d9 


^dX y© dA 
1 dr 1 d9 
VndX ~ y/QdX' 


= 9^ 


(47) 


(48) 


(49) 

(50) 


We assume the second eq.(50) is equal to /(A) and take it into the hrst equation.(49), we 
can obtain that /(A) = 1. Then these two equations will be recast into 




From f§ = 0 and = 0, we have 


(r^ + a2) {E (r2 + a^) - aL, - qQr) - aA {aE - L,) a^E^ cos^ 9 

avW 71 ' 


' + 


a 


{aLz + qQr — E {r"^ + a^)) — A(L — aE) cos^ 9 

Ay/lZ J y/9 3111^9 

When we derivative these two equations with respect to A we can obtain 

^ {E (r^ + a^) - aLz - qQr) - a{aEsin^ 9 - Lz), 

a {E (r^ + a^) — aLz — qQr) 


= 0 


^ _ f p _ 

dX [ sinH 


+ 


A 


(51) 

(52) 

(53) 

(54) 

(55) 

(56) 


Based on the above two equations, the mass of the particle plays no role but the charge of 
the particle will have physical effects. 

C. Leading order correspondence with QNMs 

Based on the WKB method, the phase of the wave function corresponds to the Hamil¬ 
ton principal function, the leading order approximation of wave equation corresponds to 
Hamilton-Jacobi equation. So the wave function has the following form: 




( 57 ) 
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By comparing to the preceding equations of QNMs, we can still find that the conserved 
particle energy E corresponds to the real part of quasi-normal frequency ojr, the z-directed 
angular momentum corresponds to the azimuthal quantum number m. When we regard 
U 0 as we can obtain Q = These is the same as what we get in Kerr back 

ground. However we cannot get the damping behaviour in leading order. In the next section 
we will discuss the next-to-leading order behaviour. 


D. Next-to-leading order correspondence with QNMs 


1. The imaginary part of quasi-normal frequency with Lyapunov exponent of radial motion 


In last section we show that the conserved quantities of particle orbit {E, Q, L^) corre¬ 
spond to the real part of parameters of a QNMs (ur, A^,m). They are the leading order 
behaviour. The damping behaviour of QNMs will appear in the next-to-leading order. In 
this section we will discuss this phenomenon. By the symmetry, we can obtain that there 
should not be any correction in 0-direction. QNMs tell us that the correction of t- direction 
will decay. Then we write the wave function as [16] 

u = AE^ = (58) 


We dehne A(t,r,9) = e '^^Ar{r)A 0 { 9 ). Then the next-to-leading order condition 29 gives 
rise to 


(dr (AdrSr) + (sinPcI,^,)) . (59) 

clT 2 \ sm9 ) 

with r the affine parameter. We should note that which is different from 

The transformation -^= can change the left hand side of this equation as: 


d In H d\nAr d In A 0 dt 
d\ ” d\ ^ '’dX' 


(60) 


The expression of || can be separated into two part: one is only depend on r and another 


on 9: 


dt - r 


(61) 


with t is the function of r and t is the function of 9. 
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rsj 


We assume that Ar is a good function near r = ro or analytic, then we have Ar{r) 
(r — ro)"" around tq with n is a integer. Then we can obtain 


7 = 1 n+- 




t 


dr/d\ 


r^VQ T 


ro {dt/d\)e 


(62) 


where 7^o = 7^"(ro). 

The physical interpretation here is that the motion in the ^-direction is independent of 
r. Then we can choose a worldline congruence whose r-value is a little bigger than ro but 
very different 6. When they return their respective initial value of 6, the change of radial r 
Ar is far less than tq. And we have that motion along a period of 6 : 

1 Ar 


= IL 


(63) 


r — ro At 

with right side will approach to a constant.As we all know the dehnition of Lyapunov expo¬ 
nent : 


A = 


lim - In 

t—>-oo t 
lAVoKO 


AX{Xo,t) 
I AXn I 


(64) 


with I AXq I represents the interval between two paths at initial time in phase space, 
I AX(Xo,t) I represents the interval between these two paths after time t. Under these 
circumstances, one path is a unstable circular orbit which the initial radial position is tq 
and will be at ro after a period At of 0 motion, the initial position of another path is a little 
bigger than ro, for example r, and will be at r -|- Ar after At. By the dehnition of Lyapunov 
exponent, we have: 

Ar 


A 


1 r + Ar — 

— In- 

At r — ro 




ro 


(65) 


1 Ar 


= 7l- 


r — ro At 

So the physical meaning of 7 ^ is the Lyapunov exponent of radial motion. The comparison 
between the circular orbit and the orbit around it tell us that the worldline will return 
their initial 6 position after a period of 6 motion. The 0-direction is Killing direction, this 
symmetry guarantee that the 0-direction of congruence will not change. Then the area of 
cross-section only depends on the radial direction: Ar(At) = r + Ar — ro- Base on the 
dehnition Ar(At) = e'^^^^Ar of the Lyapunov exponent, we can obtain 


oc e 


'TLt 


( 66 ) 
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So we have: 


A (X e 


(67) 


In the following we consider a worldline whose initial coordinates are (t, r, 6, 0) and will 
be {t + At, r + Ar, 0,0 + A0) after a period of 9. In this process, the condition f§ = 0 and 


Ig = 0 implies 


d r VW ),, 


dE 

A 

m L 


Apr* + 


A(r') 


= At; 


= 0 . 


with: 


. 6 + 


ASe = 2 ^ye{e')de' = i ^Q{e')de'. 


le. 


( 68 ) 

(69) 

(70) 


Due to the radial r is near tq and will change little, then the integral with respect to r' can be 
changed as the product of integrand and Ar. Using the Bohr-Sommerfeld condition ASg = 


2{L — \m\) and its total derivative with respect to E is zero as well as TZ{r) 
can obtain [TB] 

1 r2R"Ao 

7 = (’^+ 2 ) 7 : 


djziAL'v" 


an I an (dQ\ 

an dQ\dE) 


BS 


W, we 


(71) 


ro 


where Aq = A(ro). 

We can see that this is equivalent with the expression of imaginary part of quasi-normal 
frequency. At last we can sum up what we have done: the imaginary part of massive 
charged quasi-normal frequency of KN black hole is respond to the Lyapunov exponent of 
radial motion of massive charged particle of KN spacetime. The imaginary part of frequency 
represents the damping behaviour of wave and the Lyapunov exponent of particle represent 
the expansion of the cross section of the worldline congruence around the unstable circular 
orbits. Based on the expression (0. since IZ and Q both depend on g, the charge of the 
particle plays an important role in the motion. 


2. Angular amplitude correction and the imaginary part of Carter constant 

We can see that the Carter constant corresponds to the real part of angular eigenvalue: 

q = AL- 
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from the leading order equation, and we know that Aim is complex. We can guess that the 
Carter constant can be complex and this equation can still be hold. As we can see that 
particle charge does not appear in 0 ( 6 *), so does the angular amplitude correction Aq, so 
charge will not affect the correction. Although 0(6*) depends on the particle mass, it does 
not have essential effect on the correction. Then the meaning of complex Aim can also be 
explained by the complex Carter constant. 

From eq.( [58pp0pl| ), we can obtain the correction: 



(72) 


Ae 


All the Q and E previous are real number and in next part we will mark them as Qr and 
Er and let Q and E be complex, namely Q = Qr + iQi and E = Er + iEr As we have 
know that Ej = —7 = —w/. If we let Qi = r^) ^ind see the imaginary part as a small 

parameter and expand exp(iS'e)) which the parameters are complex, then we can obtain: 



(73) 



where the parameters in Sf are real [I 6 ] . We can see that the = Qi from the expression 


of imaginary part of angular eigenvalue. 

In summary, the correction of wave amplitude make the conserved quantity E and Q 
be complex, the imaginary part of Carter constant correspond to the next-to-leading order 


correction of angular amplitude, then the relationship Q = Aim — Tn? can hold for a complex 
Q and Aim- 


IV. EFFECTS OF CHARGE ON QNMS 


A. Unstable circular orbit and spherical orbit’s inclination angle 


We know that the radial equation : 



(74) 


is the hrst integral of motion, where : 


7^(r) = {E (r^ + a^) - aL, - qQrf - A ((L, - aEf + Q) 


(75) 
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The equation of motion is a second order differential equation, namely: 

cPr dy/lZ dr 1 dTZ 
dX^ dX d\ 2 dr 

Then the condition of unstable circular orbit means that (by [33] ) 

7^(r) = 0 = 7^'(r), 


namely 

- (A - a^) - 2a[(2M - aQ)r - Q'^]^ - Ar/ + (1 - ftl)r^ 

- 2{aQ - fllM)r^ + \a^Q^ + - /i*(a^ + Q^)]r^ 

^2c?{M - aQ)r - = 0, 


( 76 ) 


(77) 


(78) 


- (r - M)i^ + a[aQ - 2M)i - (r - M)i] + 2(1 - (79) 

— 3{aQ — jxlM)r'^ + — /i*(a^ + Q‘^)]r + a^(M — aQ) = 0. 


The conserved quantities rj = -^ and ^ ^ can be fixed by the radius r of the orbit and the 

charge parameter a = ^ of the particle. For a null curve, its parameter has a constant-scaling 
freedom A —>■ q;“^A, which induces a scaling transformation for the 4-velocity —)■ aU^. In 
fact, the equation of motion has a scaling symmetry dehned by g —)■ ag 

and —)■ aU^. The energy of the massless particle is E = g{dt, U), so the parameter a is a 
scaling-invariant ratio. Then from the particle point of view, it is the conserved qualities and 
^ that make sense. This is different from that in the QNM analysis, where the frequency, 
the counterpart of energy, is vital. 

Due to the correspondence, we know that 77(r) = 0 and 77'(r) = 0 corresponds to 
l/^(r, cu) \ro,u}R= 0 and 17'’(r, cu)' \ro,ujR= 0 in the QNM analysis, respectively. All the radii 
of spherical orbits are within ri to r 2 , and the inclination angles of these orbits will reach 
maximal values when 0 = 0. We can obtain this extreme by taking 0 = 0: 


0^(1 - IJ,^) - V - + J (a2 (1 - /i2) - 7^ - ^2)2 ^ ("X _ ^2) 


into these two equation. 
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B. Uncharged massless scalar field 

1. Redefinition of parameters 


Since the black hole is charged, the gravitational hied around black hole has been changed, 
so as the motion of particle even if the particle is uncharged. In this section, we will analyze 
this problem. In the following two sections we will see the hgure of tendency of quasi-normal 
frequency and the spherical orbit’s inclination angle with respect to the parameter of held 
as well as spacetime. In these pictures we will obtain a more direct understanding about 
the correspondence. Since there are many parameters in this problem and we are mainly 
interested in the charge ehects, we assume in the following discussions that the mass of the 
scalar held is zero. In order to simplify our discussion, we will redehne some new parameters 
hrst. 

We dehne x = Xq = ^ (representing the apex of ehective potential), y = MVLr 
(representing the real part of frequency), A = ^ (representing the angular momentum of 
KN black hole), fit = k = = ya (by the correspondence we know E = ur), p = 

and p = ^. When we consider the charged massless scalar held, the ehective potential of 
its radial equation is U’’(r) = ■^g{x,y) with 

[y{x‘^ + A^) — p\ — ycrpx]^ — (x^ — 2x + + p^) -I- ^ (1 -I- p^) y"^ — 2\py\ 

sir V) =--. 

(81) 

If we dehne 

fi = [y{x‘^ + A^) - p\- yapxf - (x^ -2x + \^ + p^) -h y (l -h p^) y‘^ - 2\py^ (82) 

and 

/2 = y[x‘^{x - 3) -h \^{x + 1) + 2p^x] + {p\ + ycrp)x - (A^ -h p^)y(Tp - p\, (83) 


then the conditions /i = 0 and /2 = 0 correspond to the requirements of QNMs, namely 
V^{r,u) |ro,a;H= 0 and W(r, ca)' \ro,LOR= 0, respectively. We can solve these two equations 
simultaneously to obtain how x and y depend on p, A, p and a. We know that ^ = 
and set h{x) = , then the imaginary part of the frequency is 


Mkli{x,y) 



2h{x)±[h{x)±g{x,y)] 

dg{x,y) 

dy 


(84) 
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2. From non-extreme to extreme 


Using numerical method, we can solve the equations([^, ([^and(84), and obtain a series 
of pictures (FigQ about how QNMs and their corresponding orbits changed when KN black 
hole approach to nearly extreme. 

We hnd that the real part of quasi-normal frequency will increase when the black hole 
charge increases, which can be found in FigjTa| Furthermore, when the black hole approaches 
nearly extreme, the relationship between y and y approaches linear when 0.788 < /i < 
l(strictly speaking, m = y = m/{l -|- 1/2) never precisely reaches ±1). What is 

important is that the imaginary part of quasi-normal frequency approach to zero at 0.788 < 
y <1 when black hole approaches nearly extreme from FigjT^ whose range is different from 
pU] . We denote y = 0.788 as yc and we call it turning point. When y > yc , the imaginary 
part of frequency will approach to zero in NEKN(nearly extreme KN spacetime). We can 
obtain that the apex of effective potential of radial equation of QNMs tends to horizon in 
NEKN, which can be found in Figj^ In Figjld[ we show the unstable spherical orbit outside 
the black hole as well as its maximum inclination angle. The horizontal axis tells us the 
radius range of spherical orbit r ^[rmin,rmax] ■ When black hole approaches to extreme we 
can obtain that many spherical orbits of different maximum inclination angle have nearly 
the same radius, r ^ 1, namely the horizon. 

Under the geometric-optics correspondence between {E, L^, Q) and (ca, m, Ai^), vanishing 
imaginary part of quasi-normal frequency corresponds to vanishing Lyapunov exponent of 
radial motion of particle in NEKN: in the aspect of QNMs, that means the mode will exist 
slightly outside the horizon for a long time and do not move away very quickly when y > yy, 
in the aspect of particle, that means a lot of particle move near the horizon with different 
maximum inclination angle and do not fall or move away after a perturbation. In addition, 
Fig(g tells us that the apex xq G [x-min, Xmax] and we can see that {xmin , Xmax) correspond 
to {rmin , ^jnax) of Same black hole parameter in FigJId] That supports the geometric-optics 
correspondence. 

It is worthy to point out that, unlike Kerr case[T6], there are inhnity ways when KN 
black hole approaches to nearly extreme. The parameters of KN(A and p) limit in a unit 
circle(A^ + < 1); so there are inhnity kinds of extreme KN black hole. We show how 

imaginary part of quasi-normal frequency changes when the parameters of KN change in 
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Real part of QN frequency{cr=0) 

y{Real part) 


Imaginary part of QN frequency{(7=0) 

z(lmaginary part) 




(a) The real part of QN frequency 


(b) The imaginary part of QN frequency 


Apex(cr=0) 

z{Apex) 



(c) The apex of effective potential 


spherical orbit's inclination angle 

cos6+ 



(d) The inclination angle of spherical orbit 


Figure 1: From non-extreme to nearly extreme. BH angular momentum A = ^ = 0.8, the dotdashed 
dark green line shows p = 0.4 , the large dashed red line shows p = 0.5, the dotted blue line shows 
p = 0.59, the thick dashed purple line shows p = 0.599, the black line shows p = 0.59999 


a unit circle with a fixed /i in Figj^ We can obtain that the imaginary part approach to 
zero when A approach to 1. Furthermore, we can obtain that not every nearly extreme KN 
black hole has the phenomenon that the imaginary part of QN frequency decay to zero. For 
example, when A = 0.4, we can require that the charge of black hole is large enough such 
that the black hole approach to nearly extreme, we obtain that the imaginary part is always 
greater than zero from Figj^ That is different from ng, In next section, we will derive the 
expression of turning point fic and obtain that fic is large than 1 in this case. 
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0.2 



Figure 2: Imaginary part of frequency 
when a — 0 and = 0.7, < 1 


Imaginary part of QN frequency{(7=0) 

z(lmaginary part) 



Figure 3: Imaginary part of frequency, 
the dotdashed dark green line shows 
p = 0.8, the dashed red line shows p = 0.9, 
the dotted blue line shows p = 0.91, the 
black line shows p = 0.9165 


3. The turning point 

In this section we require that the held is charged but massless. In NEKN, we have 
~ 1, we will solve /i = 0 and /2 = 0 under this condition, {x = l,y = is 

one of the solutions. /2 = 0 gives 


{x — l)[y{x‘^ — 2x — }? + ap) + pX] = 0. 


( 86 ) 


Then we can take y = — ^2_2x\‘i+rTp /i = 0 obtain: 

{pXY{2x - paY ={x^ -2x - X^ + apY + ^(1 + /i^)(/iA)^ 
— 2{pXY{x^ — 2x — X^ + up). 


( 86 ) 

(87) 


When considering the turning point of real part of frequency, we take x = 1 into the equation 
above and have 


^AV^ - (^A^ + p^a^ - 6pa + 6)AV^ + (A^ - pa + 1)^ = 0 

The physically relevant solution of this quartic equation is 

1 


( 88 ) 


Pc - 


V2X 


3A2 + 12 - 12pa + 2p2a2 - VE 


(89) 
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with 


S = 136 + 56A^ + - 272pa - 56X^pa + 184pV^ + UX^p^a^ - ASp'^a'^ + 4pV^ (90) 


For example, when fc = 0, we have 

ftc = i\/6(V + 4)TA^FT^iFTlS. 


(91) 


This is the turning point of QNMs of uncharged and massless scalar held in nearly extreme 


KN black hole (see H), This formula ( [9l| tells us that pc is less than 1 when A G [0.5,1], 
/ic=0.744 when A = 1 and pc = 0.788 when A = 0.8 which correspond the vanishing 
imaginary part of frequency when p > pc- But when A G [0, 0.5] ,the pc will greater than 1, 
for example pc is 1.17 > 1 when A = 0.4. That corresponds to the fact shows in Figj^ This 
agrees with the discussion in the previous section. So unlike kerr case, when the massless 
scalar helds are uncharged and the angular momentum a of nearly extreme KN black hole is 
larger than 0.5M, the imaginary part of quasi-normal frequency of scalar field will approach 
to zero with p > p^. 

Another case is that when a = 1, we can obtain that pc is limited in the range from 0 


to 0.744 by eq.(89), so when black hole turn to extreme RN black hole (A —?■ 0, p —)■ 1), the 
QNMs of charged massless held in RN have zero damping modes (See [5]). 


C. Charged massless field 

In this section we will mainly focus on the effects of field charge. 


1. Charge effects in non-extreme KN black hole 

Firstly, we want to see the influence in non-extreme KN black hole. We can observe that 
the real and imaginary parts of the frequency become greater while the apex of effective 
potential become smaller when the field charge changes from —10 to 2 from Figjla] Fig jib] 
and Figj^ When the field charge is opposite to the charge of the black hole, the electro¬ 
magnetic force is attractive like gravity, these perturbations are located at high orbits and 
possess relatively smaller energy which contain the gravitational and electromagnetic en¬ 
ergy. The electromagnetic forces change from attraction to repulsion when the field charge 
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Real part of QN frequency 

y(Real part) 



(a) The real part of QN frequency 


Apex 

z(Apex) 



(c) The apex of effective potential 


Imaginary part of QN frequency 

z{lmaginary part) 
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-- 

0.15 
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0.05 



-1.0 -0.5 0.0 0.5 1.0 


(b) The imaginary part of QN frequency 


spherical orbit's inclination angle 

COS0+ 



(d) The inclination angle of spherical orbit 


Figure 4: The influence of particle charge in non-extreme BH. The BH parameters A = a/M = 0.8, 
p = Q/M = 0.4. The dotdashed dark green line represents tr = q/E = —10, the large dashed red line shows 
a = —2, the thick dotted blue line shows cr = 0, the thick dashed purple line shows <7 = 1, and the black 

line shows cr = 2. 


changes from —10 to 2, then these perturbations need stronger gravity to balance the elec¬ 
tromagnetic force, so they are located at positions closer to the horizon and possess relatively 
higher energy. 

As we know, the particle orbits that have the maximal and minimal radii locate at the 
equator, of which the Carter constant is zero. We can rewrite the equation TZ = 0 as 

— aQr + a{a — ^) = ±g — a)A, 

A = \/— 2r + 
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(93) 





















Apex(A=0.8,p=0.4) 

z(Apex) 




Figure 5: the range of charge and apex in non-extremal KN BH 


Then we substitute this into 7^' = 0 and obtain 

- l)^^ 


a [ ± AQ — 


a ± A 


] = ±2rA - 


l)r^ 


a i y4 


(94) 


When we draw the a-r picture, we hnd that we should choose For example, see Fig,5a 
(A = a/M = 0.8, p = Q/M = 0.4). 

Fig 5a shows that the range of a has an upper bound which is saturated as r —)■ vh- 
Actually, when we set cr to be ^ in this KN black hole, we can draw the apex of the effective 
potential (Fig|5b[) in the QNM analysis and obtain that the apex will approaches the horizon 


when /i approaches 1. 

When the orbit of particle does not totally lying on the equator, or in other words, the 
Carter constant is not zero, we can not use this method and obtain the range of a easily. In 
this case, the behaviour of the worldline is much more complicated, and will be considered 
in future works. 


2. From non-extreme black hole to extreme black hole 

In this section we analyse the QNMs and their geometric correspondence when black hole 
approach to extreme. 

1. Influence of charge of black hole 

We can also obtain the phenomenon of vanishing imaginary part of frequency when 
black hole approaches to extreme with a = q/E = 1 and black hole angular momentum 
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A = a/M = 0.8, which is shown in Figj^ (The reasons why we set a = q/E = 1 and 
A = a/M = 0.8 are that we want get the charge effects and there exists zero damping modes 
in this case.) The phenomenons of real part of frequency, the apex and the inclination angle 
of spherical orbit are the same like in Figj^ One of main differents is the position of turning 
point. From the previous equation of pc we have derived, we can obtain that Hc = 0.688 


when A = a/M = 0.8 and a = q/E = 1, Fig, 6c also shows us that when 0.688 < /Xc < 1 , 
the imaginary part of frequency approach zero. The second different is that when the field 
has the same charge like black hole, the real part of quasi-normal frequency (Fig jba)) will be 
larger than that of uncharged field(FigJIa]) and the apex of effective potential(Figj^ is less 
than that of uncharged field(FigjT^. The explanation is that the electromagnetic force is 
repulsion and need more gravity to balance. So the charged perturbations will be located 
near black hole and possess higher total energy than that of uncharged perturbations. 

2. Influence of angular momentum of black hole 

When we set a fixed a = q/E = 1 and black hole charge p = Q/M = 0.8 and let KN 
black hole approach to nearly extreme (Fig |^, we can obtain the phenomenon of vanishing 
of imaginary part of quasi-normal frequency in NE KN(p = 0.8, A = 0.59999). when 0.61 < 
/i < 1 and Pc = 0.61 from FigjTbj that correspond to what we have gotten by equationjMj 

3. Disappearance of zero damping modes in NEKN 


When we use eq.(89) to calculate the turning point when A = 0.8, p = 0.6 with a = —1, 
we get Pc = 0.857 , however, when we choose A = 0.6, p = 0.8 with a = —1, we get 
Pc = 1.049 > 1. We can obtain that there is no zero damping modes in nearly extreme 
KN black hole of the last case in Figj^ Considering the variety of NEKN which is different 
from Kerr case, we show the imaginary part of quasi-normal frequency with fixed p = 1 
and (T = — 1 in Figj^ We can obtain that there are zero damping modes in a fraction of 
NEKN, but another fraction of NEKN do not have these modes when A > 0.63. In the latter 
case, the angular drag effect of the rotating black hole is not strong enough, so the gravity 
and the attractive electromagnetic force will pull the particle into black hole, resulting in 
no circular orbits or zero damping modes. So the angular momentum of black hole has 
a more important role to zero damping modes. When the angular momentum is not big 
enough under extreme case, the imaginary part of quasi-normal frequency will not approach 
to zero (See [^. 
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Real part of QN frequency{a=1) 

y(Real part) 


Imaginary part of QN frequency(cT=1) 

z{lmaginary part) 




(a) The real part of QN frequency (b) The imaginary part of QN frequency 


Apex((7=1) 

z(Apex) 



(c) The apex of effective potential 


spherical orbit's inclination angle 

COS0+ 



(d) The inclination angle of spherical orbit 


Figure 6: The influence of particle charge when KN BH approach to nearly extremel, with a = q/E = 1 
and the angular momentum of the BH ^ = -^ = 0.8. The dotdashed dark green line represents the BH 
parameters p = Q/M = 0.4, large dashed red line for p = Q/M = 0.5, thick dotted blue line for 
p = Q/M = 0.59, thick dashed purple line for p = Q/M = 0.599 and black for p = Q/M = 0.59999. 


V. CONCLUSION 

In this paper, we utilize the method of Ref.[Tn] to study the QNMs of charged massive 
scalar field in KN spacetime. We analyse the QNMs by WKB method in eikonal limit(/ S> 1) 
and get the expression of imaginary part of quasi-normal frequency first. Then we study the 
Hamilton-Jacobi equations. The charge of particle make the world line congruence be non¬ 
geodesic, but we can also find the geometric-optic correspondence under this non-geodesic 
case and obtain the corresponding relationships between the parameters of QNMs and the 
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(a) The real part of QN frequency 


(b) The imaginary part of QN frequency 


Apex((7=i) 

z(Apex) 


spherical orbit's inclination angle 

COS0+ 




(c) The apex of effective potential 


(d) The inclination angle of spherical orbit 


Figure 7: The influence of particle charge when KN BH approach to nearly extreme. (T=q/E=l, BH 
charge p = Q/M = 0.8, the dotdashed dark green line represent A = = 0.4, the large dashed red line is 

0.5, the thick dotted blue line is 0.59, the thick dashed purple line is 0.599, the black line is 0.59999 


conserved qualities of particle in leading order and next-to-leading order by identifying terms 
in the Hamilton-Jacobi equations and Teukolsky equations. We use these results to study 
the influence of charge in details. We can observe the phenomenon of vanishing imaginary 
part of quasi-normal frequency in extreme KN black hole by numerical method and then 
we obtain the expression of turning point by analytical method. The expression of turning 
point which depends on the parameters of black hole as well as charge of held gives us the 
range of zero damping modes. We also obtain a rough range of the charge of held about 
when QNMs can exist in KN black hole. 

Under the geometric-optics correspondence, the massive scalar held should correspond 
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z(lmaginary part) 
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Figure 8: Imaginary part of frequency in 


Figure 9: Imaginary part of frequency 
with /i = 1 and tr = — I 


NEKN with ct = — 1, the black line shows 


A = 0.6 and p = 0.79999, the dotted blue 
line shows A = 0.8 and p = 0.59999 

to timelike worldlines, instead of null ones, so the behaviour will be more complicated. 
Furthermore, in the case with the massive scalar held, the black hole will in general develops 
a super-radiant instability (see, e.g., [351136]), but it is not yet clear how (if possible) to see 
signals of such instability in the WKB approximation.^ This is also an interesting topic and 
should be considered in future. 
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